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Abstract 

This report provides an extensive analysis of potential wake vortex noise sources 
that might be utilized to aid in their tracking. Several possible mechanisms of aircraft 
vortex sound generation are examined on the basis of discrete vortex dynamic models 
and characteristic acoustic signatures calculated by application of vortex sound theory. It 
is shown that the most robust mechanisms result in very low frequency infrasound. An 
instability of the vortex core structure is discussed and shown to be a possible mechanism 
for generating higher frequency sound bordering the audible frequency range. However, 
the frequencies produced are still low and cannot explain the reasonably high-pitched 
sound that has occasionally been observed experimentally. Since the robust mechanisms 
appear to generate only very low frequency sound, infrasonic tracking of the vortices may 
be warranted. 
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Introduction 


Wake vortices present a potential hazard to aircraft in the vicinity of airports. For 
this reason, aircraft operations are scheduled so as to permit the necessary time for the 
wake vortices of preceding aircraft to either disperse or completely move out of a 
protective zone based on a worst case scenario. This results in very inefficient use of the 
existing runway capacities around many airports. If the position of the vortices were 
known at all times, aircraft operations could then be carried out much more safely and 
efficiently. One scheme for actively monitoring such vortices is to passively track the 
acoustic emissions that they produce. This possibility is currently under examination 
both in terms of the fundamental wake acoustics phenomenology and the operational 
concept involving such a scheme by the U.S. Government. 

In previous acoustic studies within the program, four candidate source 
mechanisms had been identified to be appropriate for further concentrated examinations: 

1 . Core vibrations excited by initial conditions or instabilities 

2. Unsteady transport of the vortex, including possible interactions with 
secondary vortices shed from the atmospheric boundary layer 

3. Vortex core motion excited by intermittent turbulent structures 

4. Unsteady advection of turbulence around the vortex core 

It was also stated that the goal of the studies was to determine whether a unique 
signature (one that can be identified in the presence of noise) of a consistent dynamic 
characteristic (one that is present for wake vortices during the time of interest for 
detecting them) exists. The purpose of this report is to address the aforementioned 
question. This objective is accomplished by applying classical aeroacoustic theory to 
engineering vortex dynamic models of these wake phenomena, in order to determine 
whether a reliable acoustic signal can be identified for tracking purposes. 


Vortex Dynamic Models 

In order to investigate the sound generation by aircraft wake vortices, engineering 
models of the phenomena based on the Lagrangian discrete vortex method are developed. 
Since these phenomena are inherently vortical in nature, vortex method is therefore 
particularly well suited and very efficient in describing the fundamental fluid mechanics. 
It is worth noting that engineering models of this type have demonstrated their utility in 
studies such as NASA’s AVOSS program [1], In addition, these fluid dynamic results 
can be conveniently integrated with the theory of vortex sound to identify the 
characteristic acoustic signature associated with the phenomenon. 
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a. Wake Vortices Interacting with the Ground Plane 

The logical first mechanism to explore would be the sound generated by the 
motion of the vortex pair. After aircraft vortices are fully rolled up, the port and 
starboard vortices mutually induce each other to move downward until they feel the 
presence of the ground. As they approach the ground, the increased pressure produced by 
their interaction with the ground causes the vortices to separate and propagate outward 
from the aircraft’s path. As a model of this phenomenon, consider an initial 
configuration of two wake vortices of strength r and opposite senses of rotation 
generated by an aircraft of span 5 at an altitude h as shown in Figure 1 : 


-r r 



It may be of interest to note that experimental evidence has revealed that at a 
given location along the flight path, the fundamental vortex dynamics can be described as 
a two-dimensional problem [2], as modeled in Figure 1. The initial separation, b Q , of the 

vortices has been found to be given by b Q = (n / 4) s . The strength of these vortices may 

also be estimated since the lift per unit span of the aircraft is given by / = p U y , where 
p is the density of the air, U is the speed of the aircraft and y is the circulation per unit 
length in the wake of the wing. For equilibrium or quasi-equilibrium flights, which are 
representative of aircraft cruise and landing configurations, the total lift will equal the 
weight of the aircraft, W . Thus, 


4 W_ 
n p U s 


( 1 ) 


where the factor 4 / n stems from the assumed elliptic lift distribution on the wing. 

Consider a coordinate system such that one vortex lies at (x,y) and the other at 


5 



(- x,y ) as shown in Figure 2: 


y 



Figure 2: Wake Vortex Coordinate System. 


To determine the vortex trajectories, note the initial conditions 


x(0 ) = (tt / 8) 5 

y( 0 ) = h 


for the vortex in the first quadrant. The presence of the ground plane can be handled by 
inserting image vortices as shown in Figure 3: 


y 
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which satisfy the boundary condition that the normal component of the induced velocity 
field is zero on the ground plane. Then, applying the Biot-Savart Law [3], the velocity 
components induced on the vortex core in the first quadrant are: 


dx _ r x 2 
dt Any (x 2 + y 2 ) 

dy_ _ -r y 2 
dt 4 K x (x 2 + y 2 ) 


This coupled set of nonlinear, first order ordinary differential equations can be 
simply numerically integrated starting from the initial conditions given by Equation (2) to 
determine the vortex locations as a function of time. The path of the vortex in the second 
quadrant is obtained by symmetry. Note that the case of a crosswind can be handled in 
exactly the same way, by assuming the lateral transport of the vortices advect with the 
prescribed crosswind profile. The only difference is that the vortices do not remain in 
mirror image positions and thus a set of four equations must be solved to track both 
vortices. 

As an example under zero crosswind, the case of a B747 with wingspan of 59.7 
m, altitude of 150 m, and a constant circulation strength of 595.54 m 2 /sec. has been 
considered. The time step, At = 0.5 sec, was employed in conjunction with a second- 
order numerical integration method. Figure 4 shows the trajectory of the vortex core in 
the first quadrant in a time frame of 150 seconds. 
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The path of the other vortex is the mirror image about the y-axis. The trajectories 
approach asymptotically to an altitude of 23.4 meters, which represents a height of 
b a / 2. When the wake is at the aforementioned height, it is considered to be fully in 
ground effect. Figure 5 displays the vortex positions, x and y, while Figure 6 displays the 
vortex velocity components, u = dx/dt and v-dy/dt, for the vortex in the first quadrant 
as a function of time after the vortices are shed from the aircraft. As time goes to infinity, 




Figure 6: Vortex Core Transport Velocity Components. 


the v velocity (i.e., vertical component), which initially was -1.97 m/sec goes to zero and 
the u velocity (i.e., horizontal component) asymptotically approaches to a constant value 
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of 1.99 m/sec., which represents the ground-induced self-propelling vortex velocity. 
Note that these u and v limiting values take place at the wake altitude of b 0 /2, and the 

particular numerical example showed that the process took 70 seconds for the wake to be 
fully in ground effect. The importance of this time frame will become evident in the 
subsequent acoustic analysis. Finally, Figure 7 displays the acceleration components 
du / dt and dv / dt . Both of these acceleration components approach asymptotically to 
zero as time progresses. These acceleration components are important to the vortex noise 
generation as will be seen later. 



Figure 7: Vortex Core Acceleration Components. 


Note that the example has been computed using constant circulation values for 
both portside and starboard vortices. The process of vortex decay, which would alter the 
acceleration characteristics, can be prescribed from auxiliary equations. Additional 
fidelity of the vortex dynamics may include the effects of altitude wind shear, 
atmospheric boundary layer profile and vortex induced secondary vorticity, whose net 
effects can all be modeled as discrete vortices. The flow field is then computed using the 
generalized form of the Biot-Savart Law, by examining the induced velocities from all of 
physical and the associated image vortices. The associated modeling and numerical 
issues are beyond the scope of the present study. The aforementioned description is 
provided for the purpose of emphasizing that higher order wake modeling can be 
accomplished within the present theoretical framework. 

b. Two Unequal Magnitude Vortices in Free Space 

The problem of two unequal magnitude vortices of the same sense of rotation in 
free space was also considered as a model representing one of the wake vortices 
interacting with an atmospheric eddy, a wind shear, or other airframe vortices such as that 
of a flap. Consider the coordinate system shown in Figure 8: 
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Figure 8: Unequal Vortices Coordinate System. 


i.e., a vortex of strength r x at the position (x, ,y x ) and a vortex of strength r 2 at the 
position (x 2 ,y 2 ). Again applying the Biot-Savart Law, the equations governing the 
system are 


dx x 

dt 

d}\ 

dt 

dx 2 

dt 

dy 2 

dt 


£2 

2 nr 1 

-_h 

2 nr 1 

zh 

2 nr 1 
2 nr 1 


O2-L1) 

(x 2 — Xj ) 

(V2-L1) 

(x 2 -Xi) 


where r 2 = (x 2 - x, ) 2 + (y 2 - y t ) 2 can be shown to be a constant equal to its value at the 
initial condition. This set of equations may readily be solved by Laplace transformation 
yielding the solution 


x, ( t ) = X[ (0) + X x (t) cos 6 - 7j ( t ) sin 0 
v, (0 = y, (0) + X l ( t ) sin 6 + Y l (t) cos 6 
x 2 (t) = Xj (0) + X 2 {t) cos 0 - Y 2 (t) sin 6 
y 2 it) - y x (0) + X 2 (t) sin 6 + Y 2 ( t ) cos 0 
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where cos# = (x 2 (0) - x, (0 ))/ r , sin 0 = {y 2 ^)-yi^))lr, and 


A + B A+B 
Y [ (t) = ^-sm(A + B) t 


X 2 {t) = 


B r 


■ + ■ 


A r 


A+B A+B 


cos(t4 + 5) t 


Y 2 {t) = ^-sm{A + B) t 

A + B 


where A = r x /(2 n r 2 ), B = T 2 1(2 n r 2 ), and the values at time = 0 are the initial 
conditions. Note that the solution is periodic with frequency co = A + B or 
/ = (A + B)/(2 n) as the vortices merely rotate about one another periodically. Since 
the propagation of acoustic waves is a linear phenomenon, this same frequency will 
appear in the acoustic field. Some idea of the frequency of this oscillation may be 
obtained by again considering the case of a B747 where /j =595.54 m 2 /sec. Supposing 
for example that _T 2 = 0.1 • /j and r = 10 m yields / = 0.16 Hz, a frequency much below 
human hearing (30 Hz). Even if r = 1 m, which would realistically be much too close, 
the frequency would be only 16 Hz. Thus, this mechanism produces only infrasonic 
signals and could not be responsible for any periodic signal that might be heard. 

As was seen before, the presence of a ground plane can be added into this model 
quite simply using image vortices. However, then the motion can no longer be solved 
analytically. The solution is known [4] from previous numerical calculations to exhibit 
chaotic behavior. Such behavior is not a numerical artifact but is produced by extreme 
sensitivity to initial conditions. For some sets of initial conditions, the two vortices will 
continue to oscillate periodically about one another as they are convected by their images 
in the plane as was seen in the previous analysis of wake vortices in the presence of the 
ground. For other initial conditions, the two vortices separate as each is convected at 
different speeds by their image vortices. 


Applications of Vortex Sound Theory 

The theory of vortex sound generation was first developed by Powell [5] who 
showed that the primary source of sound generation in low Mach number flows was due 
to the presence of vorticity, in fact, acceleration of vorticity is the major source. The 
acoustic pressure, p , was shown to satisfy the wave equation 
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( 5 ) 


1 d 2 p 
c 2 dt 2 


-V 2 p = p 0 V -(fflxv) 


Here c and p o are the ambient speed of sound and density respectively and co and v are 

the vorticity and velocity vectors. Integrating Equation (5) and performing a multipole 
expansion, which amounts to neglecting retarded time differences across the source 
distribution, it can be shown [6] that the solution to this equation in free space is given by 


p(X,t) 


Po x i x J 


( 


4 n c 


X 


d 2 
dt 2 


J 




Lj dy 


V 


( 6 ) 


Here X = (x l ,x 2 ,x 3 ) is the observer position, L = coxv is the Coriolis acceleration 
vector, V is the volume containing the flow, and the asterisk indicates evaluation at the 


retarded time 


t- 


X 


/ c . 


a. Wake Vortices in Presence of Ground Plane 


When the flow takes place in the presence of a surface, one ordinarily expects an 
additional surface integral [7] to appear in the solution given by Equation (6). However, 
for the special case of a plane surface, Powell [8] showed that the surface integral is 
negligible if the image vortices are treated as real and the volume integral taken over both 
the real and image spaces. Thus, consider the geometry shown in Figure 3 for the case of 
wake vortices interacting with a ground plane. The vortices have velocities only in the 

y x and y 2 (or x and y) directions, i.e., for the vortex in the first quadrant v = (u,v, 0) 
where u and v are shown on Figure 6. Further, they only have a vorticity component in 
the y 3 direction, i.e., again for the vortex in the first quadrant 0) - (0,0, rS(y - Tv)) 
where c)(«)is the two dimensional Dirac delta function 8(x- x v )8(y - y v ) and 
y v - (x v ,y v ) is the vortex position in the x — y plane. Thus, the Coriolis acceleration 
only has components in the y x and y 2 directions, i.e., L - (~rv8(y - y v ),ru8(y - y v ), 0) 
for the vortex in the first quadrant. Thus, jy ( L 3 dy - 0 and, assuming that the coordinate 

V 


system is placed such that the length l of the vortices extends over -l/2<y 3 <1/2, 



= 0. This choice for the origin of the coordinate system not only simplifies the 
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integrations, but also minimizes the retarded time differences over the length of the 
source distribution. Carrying out the other 4 integrals over all four (real and image) 
vortices, paying attention to the positions and signs of both the circulations and velocity 
components yields 



= -4 xfvl 



V V 


J y 2 L 2 dy - Ayful 


v 


where x and y are the vortex positions shown in Figure 5. Thus, applying these 
relations in Equation (6) yields the expression for the acoustic pressure radiated by the 
vortices 


p(X,t + \X\/c) = -P^ 
K C 2 X 


at at 


( 7 ) 


Now, suppose that the observer is on the ground at a distance d from the aircraft 
flight path and at the center of the vortices’ length, i.e., X - (d, 0,0) . Then the 

expression for the observed acoustic pressure, given by Equation (7), becomes: 


p(X, t + d / c) = 


PoSl^l 

n c 2 d dt 2 


(xv) = 


Po£ j 

K c 2 d 


. dv d 2 v du 

2 U b X T- + v — 

dt dt dt 


In order for the use of this expression to be precisely valid, the observer should be 
in the farfield of the source, i.e., the distance d should be large enough that both d > / 
and d > X , where X is the wavelength of the sound produced. Both of these conditions 
are often violated in wake vortex noise monitoring. However, all the violation means is 
that higher order sources may come into play and that the amplitude of the acoustic wave 
may not fall off like 1 / r . The basic signal characteristics, which are very low in 
frequency as discussed below, will be preserved. 

The normalized acoustic pressure at the observer time t + d / c is given by 
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( 8 ) 


n c 2 d 

p. r i 


. d 2 . . . dv d 2 v 

p{X, t + d/c )- — t (xv)-2u hx — — 

dt dt dt 


du 


+ v — 

dt 


As an example, the normalized acoustic pressure radiated in the case of a B747, as 
discussed above, is plotted in Figure 9 as a function of source time t rather than observer 
time t + d/c for an observer at a distance of d = 50 in from the ground track of the 
aircraft and an initial vortex altitude of h = 150 m. The signal has a zero asymptote as 
time continues and the acceleration of the vortex diminishes. Note that the acoustic 
signal resembles one cycle of a periodic signal with a period near 45 seconds that would 
correspond to a frequency near 0.02 Hz, again inffasonic. 



b. Elliptic Vorticity Distribution 

As an example of a distributed vortex, consider a uniform distribution of vorticity 
Q over an elliptical region in free space as shown in Figure 10: 
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Wake vortices of this shape have been confirmed via smoke visualization flight tests and 
this appearance might be due to the influence of the fuselage or wing on the rolling up of 
wakes. The origin of coordinates is again taken at the center of the vortice’s length with 
the coordinate y 3 in the axial direction. Thus the vorticity vector is given by co = Q k . 

This vorticity distribution in the infinite length, two dimensional case has been studied by 
several research workers, including Howe [9], who have shown that the configuration is 
stable and rotates with angular velocity dO / dt = Q / 4 to the order of s , where s is the 
eccentricity of the ellipse. Due to the rotation, the boundary of the vorticity distribution 
at any time t is given by r - R-a[ 1 + e cos(20 -(Q / 2)t)\. The velocity vector lies in 

the plane, i.e., v = u i + v j , where 


u — 


v — 


— £2 v [ 

sin 6 + e sin(0 


2 

Q r 


Q 

cos 6 - e cos(0 - — t) 


Thus, the Coriolis acceleration vector becomes L - -Q v i + Q u j , i.e., L - (L X ,L 2 , 0) . 
The sound radiated by this vortex motion may be calculated from Equation (6). 


15 



Since L 3 = 0 , jy i L 3 dy = 0. Further, due to symmetry, jy 3 Ljdy = 0. Thus, the 

V V 

necessary integrals, correct to 0(e) , are: 



V 



L£2 2 a 4 x., Q. . 

(1 + e cos — t) 

8 2 


LQ 2 a A n .. Q x 

(1-ecos — n 

8 2 


and 



v v 


LQ 2 a 4 K e . Q 

sin — t 

8 2 


where L is the length of the vortex. The observer will again be taken to be in a plane 
through the origin which lies at the center of the vortex, i.e., X - (x, , x 2 ,0) , as shown in 

Figure 10. Utilizing the above integrals in Equation (6) yields an expression for the 
sound radiated to the observer position 


p( Y, l) = z£Mddl 

128 c 2 X 


Q 

T 


(Xj 2 -x^cos — 1 + 2xjX 2 sin — t 


Q 

T 


where, again, the asterisk indicates evaluation at the retarded time t — 


X 


/ c . This 


expression may be simplified somewhat by introducing the polar coordinates 


X 


cosy/ and x 2 = Ixlsini// yielding 


p(X, t + 


X 


. , - p LSI a e 

/ c) = — — cosl 


128 c 2 


X 


( Q 

2y/-—t 


( 9 ) 


Thus, it can be seen that the sound radiated by the rotation of the elliptical 
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vorticity distribution is periodic with frequency / = Q / An . Note that the amplitude of 
the acoustic pressure is proportional to the ellipse eccentricity s . 

A concept of the importance of this mechanism may be obtained by considering 
modeling a typical wake vortex by this elliptical vorticity distribution. The area of the 
ellipse is A = n a 2 , again to the order of s , and the circulation around it is given by 
r - Q A . Thus, for modeling this situation, it is necessary to take 



r 

— 2 

n a 


yielding a frequency that is correct to the order s . 


f = 


£ 2 _ 

4k 


r 

4k 2 a 1 


As an example, consider a B747 with r = 595.54 m 2 /sec. and a = 1.19 m, then 
/ = 10.76 Hz. Thus, it can be seen that this elliptical vortex rotation mechanism also 

produces very low frequencies for typical wake vortices. 


Hydrodynamic Pressure Field 

Due to the fact that the amplitude and frequency of the sound calculated above 
were very low, it is of interest to compare this sound field with the hydrodynamic 
pressure field produced by the vortex. A vortex produces its own pressure field which is 
carried along with it as it moves. This hydrodynamic pressure field is not acoustic, i.e., it 
does not propagate as a wave, but will cause a significant pressure fluctuation as it passes 
over an observer location. For example, it can be shown that the pressure field produced 
by a line vortex in free space is given by 
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pir) = Po~ 


r 2 


8 K 2 r 


2 


where r is the radial distance from the vortex core center and p 0 is the ambient pressure 
(i.e., stagnation pressure of the atmosphere). Note that the pressure near the vortex is 
lower than ambient and returns to ambient like Mr 2 . 

For the phenomenon of wake vortices near a surface, for simplicity only one 
vortex and its image are considered as shown in Figure 11. In addition, zero crosswind is 
assumed. 


y 



Figure 11: Single Wake Vortex and Its Image. 


This model is based on the consideration that over times of interest the vortex will 
be more influenced by its image than by the wake vortex from the other end of the wing, 
as its induced velocity becomes negligible. In this case, the vortex is convected by its 
image at a constant convection velocity, U c - T /(4 n y v ) where y v is the constant 
altitude of the vortex. Recall the example for the B747 shown in Figure 6, where y v will 
become b 0 /2 and U c approaches to 1.99 m/sec. The observer will be assumed to be on 
the surface at a distance r from the vortex where the induced velocity components are 
u = (r y v )l(n r 2 )and v = 0 in the x and y directions respectively. The pressure field 
on the surface may then be shown to satisfy the partial differential equation 


d_ 

dx 



= 0 
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Integrating this equation yields an expression for the deviation of the gauge pressure from 
ambient on the surface 


Ap(t) = p-p o = p o uU c - 


( 10 ) 


Equation (10) has been calculated for the case of a B747, again for an observer at a 
distance of d = 50 m from the ground track of the aircraft, with an initial vortex altitude 

of h = 150 m. and at the center of the vortice’s length /, i.e., X = (t/,0,0) . The time 
history of the hydrodynamic gauge pressure field produced by this calculation is shown 
in Figure 12, nondimensionalized by the same parameters as were used for the acoustic 
pressure. The vortex length has been taken as l = 1000 m in the nondimensionalization. 



Figure 12: Hydrodynamic Pressure Time History. 


Note that the gauge pressure fluctuation is initially negative as the observer feels the low 
pressure near the core of the vortex. Then, as the vortex passes overhead of the observer 
position, the pressure fluctuation rises dramatically to an overpressure that then dies off 
to zero as the vortex is convected away from the observer location. In comparing this 
hydrodynamic pressure time history with that of the acoustic time history for the same 
aircraft shown in Figure 9, it can be seen that the hydrodynamic pressure field is several 
times larger than the acoustic field at this observer location. Further, the hydrodynamic 
pressure fluctuations occur over nearly the same time interval since the acoustic delay 
time, die, is only 0.15 seconds. Of course, the acoustic field only falls off with distance 
like 1/r while the hydrodynamic field falls off like Mr 1 . Thus, the acoustic field will 
eventually become dominant. However, for observer positions such as being used for 
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monitoring wake vortices near the ground, the hydrodynamic pressure field will 
dominate. 

The dominance of the hydrodynamic pressure field and characteristic signature 
shown in Figure 12 suggests that measurements of the hydrodynamic pressure field might 
well be employed to track the vortices. As can be seen, the fundamental frequency of the 
hydrodynamic pressure field is near 0.01 Hz, similar to that of the acoustic pressure field. 
Pressure fluctuations with frequencies of this order are readily measured by means of a 
microbarograph [10] such as are being used to monitor nuclear testing throughout the 
world. Such pressure transducers are fairly large, on the order of a foot in diameter, and 
typically have a frequency response that is flat over the interval 0.02-5 Hz. Their 
response is only 3 dB down at the frequency of interest here, / = 0.01 Hz. 


Instabilities 

Another mechanism that has been suggested as a possible source of wake vortex 
noise is the presence of instabilities within the vortex itself. If the vortex were purely 
two-dimensional, Michalke and Timme [11] have shown that it would be stable by the 
Rayleigh criterion. However, aircraft wake vortices actually contain a significant axial 
velocity component, which is explained by first by Batchelor [12] as arising from the 
competition of circulation-induced pressure gradient and entrained viscous fluid shed by 
the wing boundary layer into the core. The shear in the axial velocity combined with the 
vorticity-induced swirl velocity can cause the wake vortex to be unstable to certain types 
of disturbances as proved by Lessen, Singh and Pallet [13]. They found the vortex to be 
unstable to non-axisymmetric disturbances with negative azimuthal wavenumbers. The 
theory of Lessen, Singh and Pallet has been validated by Khorrami and Singer [14] using 
measurements of the wake vortex produced by a part-span flap. Although they found the 
instabilities to grow for all negative azimuthal wavenumbers, the data indicated that those 
with the peak amplification rates, i.e. those in the frequency range 


2 < 


(Q r_o 

U 


<4 


( 11 ) 


Here U is the mean velocity in the streamwise direction, r is the core radius of the 
vortex and co is the frequency in radians per second appeared to dominate. 

This relation can be applied to the case of aircraft wake vortices in order to 
estimate the frequencies that might be produced by this instability mechanism. The best 
data available on wake vortex core radii indicates that the radius may be approximated as 
r 0 - 0.02 .s' , where s is, again, the span of the aircraft [15], Taking U = 66 m/sec as a 

typical landing speed, Equation (11) yields the following estimates for aircraft in the 
commercial fleet: 
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Table: Predicted Core Radii and Instability Frequencies 


Aircraft Type 

Core Radius 

Frequency Range 

B747 

1.19m 

17.6-35.2 Hz 

B757 

0.76 m 

27.6-55.2 Hz 

B767 

0.95 m 

22.1-44.2 Hz 

MD80 

0.66 m 

31.9-63.8 Hz 

A320 

0.68 m 

31.0-62.0 Hz 


Since the sound radiation produced by such instabilities satisfies a linear wave equation, 
the sound radiation will occur at the same frequency as the instability. Note that these 
predicted frequencies, while mostly audible, are still on the low end of the normal range 
of human hearing. 

The instability mechanism discussed above occurs in incompressible flow. There 
is also the possibility of a compressible instability since the vortices are formed at altitude 
where the density is different than that nearer the ground. Since, in the absence of 
viscous action, vortical fluid tends to remain vortical and vice versa, the vortices would 
carry the density acquired at altitude with them as they descend allowing buoyancy 
effects to come into play. Further, it is known that a Kelvin-Helmholtz instability can 
occur between fluids of differing density and velocity. Such potential instability 
mechanisms have not been examined as yet. 


Conclusion 

This report has presented a technical review of the classical theory of vortex noise 
generation as applied to wake vortex noise generation in support of the wake acoustics 
research program. Engineering vortex dynamic models were developed for several 
phenomena that have been considered as possible candidate mechanisms for producing a 
reliable acoustic signature for tracking purposes. All of the mechanisms considered 
resulted in low frequency sound generation and suggest that infrasonic tracking of the 
vortices might be profitably employed. Higher frequency sources would require fluid 
mechanical mechanisms that oscillate a high frequency, possibly suggesting small scale 
turbulent phenomena. Such phenomena will be the subject of further research. 
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